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The dissipation associated with nonequilibrium flow processes is re- 
flected by the formation of strange attractor distributions in phase 
space. The information dimension of these attractors is less than that 
of the equilibrium phase space, corresponding to the extreme rarity of 
nonequilibrium states. Here we take advantage of a simple model for 
heat conduction to demonstrate that the nonequilibrium dimension- 
ality loss can definitely exceed the number of phase-space dimensions 
required to thermostat an otherwise Hamiltonian system. 

PACS: 05.20, 5.45.D, 05.70.L 



Nonequilibrium molecular dynamics has been used to establish a close link between 
microscopic dynamical phase-space instabilities and the macroscopic irreversible 
dissipation associated with the Second Law of Thermodynamics [1,2]. Both non- 
Hamiltonian and Hamiltonian methods have been used [3] . The simplest such con- 
nection between microscopic dynamics and macroscopic dissipation results when 
one or more Nose-Hoover thermostats [4] are used to control nonequilibrium steady 
states. In the absence of nonequilibrium fluxes and with sufficient phase-space mix- 
ing, these thermostats generate Gibbs' and Boltzmann's canonical distribution. In 
cases which include nonequilibrium driving the instantaneous external entropy pro- 
duction rate (due to heat transfer with the thermostats) is proportional to the sum 
of the instantaneous Lyapunov exponents [5]: 

S/k = * = ~®/ ® ■ 

The summed-up instantaneous exponents give the comoving rate of change (fol- 
lowing the motion) of an infinitesimal phase-space hypervolume <g>. The sum of 
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the time-averaged Lyapunov exponents, {(A)} is necessarily negative, reflecting the 
irreversible dissipation described by the Second Law of Thermodynamics. 

Despite the time reversibility of the underlying equations of motion, a symmetry 
breaking due to the enhanced stability of trajectories proceeding "forward" in time 
relative to their reversed twins, guarantees irreversibility and the formation of a 
phase-space strange attractor [5]. The dimensionality of such a strange attractor 
can be measured by computing the smallest number of time-averaged exponents 
(beginning with the largest one, (Ai)) with a negative sum. This "Kaplan- Yorkc" 
dimension is identical to the "information dimension" in cases of physical interest 
[6]. 

Because Liouville's Theorem [7] establishes that the "extension in phase," the 
hypervolume of the comoving element ®, cannot change in a motion governed 
by Hamiltonian dynamics, there has been some reluctance to accept this strange- 
attractor explanation of irreversibility. Very recently Ramshaw has provided a par- 
ticularly clear analysis of the generalized Liouville's Theorem required to describe 
noncquilibrium dynamics with time-reversible thermostats [8]. Refs. 9 and 10 pro- 
vided a useful computational model and a means to estimate the dimensionality 
loss, AD, the difference between the strange attractor's information dimension and 
that of the phase space in which it is embedded. It was shown that an accurate 
estimate for AD can be made, under certain conditions, which requires only a single 
exponent, not the whole spectrum. We use the whole spectrum here because we 
want to make the demonstration of dimensionality loss as convincing as possible. 

In the "</> 4 " model an otherwise harmonic nearest-neighbor lattice, with the 
Hooke's Law pair potential for Particles i and j, 

0(r) = i(r- l) 2 ; r = \r i -r j \ >0, 

has each particle tethered to its lattice site with a quartic potential, jSr 4 . The 
quartic tethers prevent momentum conservation, so that Fourier heat conduction 
can be observed, and also can provide chaos, with one or more positive Lyapunov 
exponents. To model a nonequilibrium heat-conducting state two or more of the 
particles arc thcrmostatcd using feedback forces {~Cp} which are linear or cubic in 
the friction coefficient £ and the momentum p. Here we choose the simplest case, 
linear in both variables, the Nose-Hoover thermostat. In one space dimension the 
thermostated equations of motion are: 

p = F-(p ;C=[(p 2 / m kT)-l}/r 2 , 

where r is the characteristic response time of the thermostat force, — Qp. In two 
dimensions, where p 2 is a sum of x and y components, T is replaced by 2T. The 
target temperature T for the thermostat is necessarily achieved whenever a station- 
ary solution exists, such that the long-time-averaged value of £ vanishes. Because 
the long time average of (( — (d/dt)^( 2 likewise vanishes, the heat transferred by 
any thermostated momentum can be directly related to the temperature: 

((p 2 /m) = (C)fcT , 

again explicitly writing only the one-dimensional case. 

The I6-body system shown in Fig. 1 has a cold particle at the lower left corner, 
a hot particle at the upper right corner, along with 14 Newtonian particles able to 
transmit heat from the hot particle to the cold one. This system (as well as many 
others) can exhibit dimensionality losses exceeding the 10 phase-space coordinates 
{x, y,p x ,p y , (} required to describe the two thermostated particles. In the case 
shown in the Figure, with cold temperature 0.001 and hot temperature 0.009, and 
a Nose- Hoover relaxation time r = 1.4, the dimensionality loss is 12.5. This means 
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FIG. 1. Stationary state particle trajectories for a 16-particle system with dimension- 
ality loss of 12.5 induced by two Nose-Hoover thermostated particles at temperatures of 
0.001 and 0.009. 

that the sum of the 53 largest Lyapunov exponents, plus half the 54th, is zero. See 
Fig. 2 for the spectrum of exponents. The strange attractor has an information di- 
mension of 53.5, embedded in a 66 dimensional phase space, of which 56 dimensions 
represent the purely Hamiltonian particles. 

This result establishes very clearly that the dissipation and overall contraction 
occurring in the thermostated part of the system can cause a loss in the Hamiltonian 
region too. Although we had been able to simulate systems with dimensionality 
losses barely exceeding the additional boundary phase-space coordinates [11], the 
present results are much more clearcut. Otherwise one might imagine that the 
thermostat regions simply provide a time-dependent force, which for Hamiltonian 
dynamics can provide no dimensionality loss. The feedback linking the forces to the 
phase-space coordinates, not just to the time, makes the loss possible. 

This situation is analogous to the dynamics of a one-dimensional damped har- 
monic oscillator with a fixed friction coefficient z, for which the logarithmic phase- 
space contraction rate, 



The damped oscillator motion nevertheless collapses (<8> — > 0) to the fixed point at 
the origin, (q = 0,p = 0) because the motion rotates the comoving phase volume. 
Phase-space rotation [12] in many-body systems increases rapidly with system size 
and is necessarily the mechanism through which the volume collapse can spread to 
directions without direct contraction. 

In summary, thermostated heat flow, in a tethered harmonic lattice, shows conclu- 
sively that the dimensionality loss associated with irreversible processes can exceed 



®/<8> = (dq/dq) + {dp/ dp) , 



occurs in the momentum direction only: 



q = p\ p=-q- zp 



(dq/dq) — ; (dp/ dp) = —z . 
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FIG. 2. Lyapunov spectrum for the system shown in Fig. 1 There are 24 positive 
exponents, a single vanishing exponent (emphasized here), and 41 negative exponents 
describing the motion in the 66-dimensional space. In the symmetric equilibrium situation 
(with the hot and cold temperatures equal) the exponents occur in 33 pairs, {±A}, with 
a single pair of vanishing exponents. In the nonequilibrium case there is an overall shift 
toward more negative values, S/k = — ^2 A > 0. 

the additional variables required to specify the thermal boundaries driving the sys- 
tem from equilibrium. This shows that the dimensionality loss is a real feature of 
nonequilibrium systems, not an artifact of a particular choice of thermostat. We 
believe that this finding clarifies a point which has been strenuously debated over 
the past decade [13]. There can be no doubt that advances in parallel computing 
will soon elucidate the precise way in which the large-system limit is obtained. 

The approach we follow here connects irreversible entropy production to mul- 
tifractal phase-space structures through time-reversible Nose-Hoover thermostats. 
It is worth mentioning that there are now alternative connections, including some 
based on purely Hamiltonian mechanics [3]. 
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